Implementing optimal control pulse shaping for improved single-qubit gates 
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We employ pulse shaping to abate single-qubit gate errors arising from the weak anharmonicity 
of transmon superconducting qubits. By applying shaped pulses to both quadratures of rotation, a 
phase error induced by the presence of higher levels is corrected. Using a derivative of the control 
on the quadrature channel, we are able to remove the effect of the anharmonic levels for multiple 
qubits coupled to a microwave resonator. Randomized benchmarking is used to quantify the average 
error per gate, achieving a minimum of 0.007 ± 0.005 using 4 ns- wide pulse. 

PACS numbers: 03.67.Ac, 42.50.Pq, 85.25.-j 



The successful realization of quantum information pro- 
cessing hinges upon the ability to perform high-fidelity 
control (gates) of quantum bits, or qubits. A value of 
10 -3 error per gate (EPG) is typically quoted as the nec- 
essary threshold for fault-tolerant quantum computation 
PQ . For any two- level system, the optimal achievable gate 
performance is set by the ratio of gate time to coherence 
time. However, quantum information processing systems 
consist of multiple coupled qubits. Often these qubits are 
not truly two-level systems, but rather quantum objects 
with a rich level structure. Thus, one important chal- 
lenge is to achieve optimized single-qubit gates in a large 
Hilbert space. 

Optimal control theory has been previously employed 
in nuclear magnetic resonance, non-linear optics, and 
trapped ions to combat specific errors such as always- 
on qubit couplings and spatial inhomogeneities [2J. In 
these systems, even with optimized qubit control, lim- 
its to qubit gates are often due to systematic errors 
rather than decoherence. Superconducting qubit control, 
however, has primarily been limited by coherence times. 
With recent progress such as the demonstration of high- 
fidelity single-qubit gates [3J 2|, two-qubit gates [ME], 
entanglement [5141 1| and simple quantum algorithms [7], 
the superconducting qubit architecture is growing into 
a more complex quantum information testbed, placing 
the level of qubit control under increased scrutiny. One 
approach towards improving single-qubit gates is to de- 
crease the total gate time. However, in multi-level qubits 
such as the transmon [T5] or phase qubit [3J, the weak 
anharmonicity sets a lower limit on the gate time. Fur- 
thermore, superconducting qubit coupling schemes such 
as circuit quantum electrodynamics (QED), in which a 
transmission- line cavity couples multiple qubits [131 114] , 
can make single-qubit control more difficult as a result of 
many extra levels in the Hilbert space. 

In this Letter, we implement simple optimal control 
techniques to improve single-qubit gates in a circuit QED 
device with two superconducting transmons. The pulse- 
shaping protocols we investigate are guided by the recent 



theoretical exploration of derivative removal via adia- 
batic gate (DRAG) in Ref. [15] . We demonstrate the im- 
provement of single-qubit gates on two separate qubits 
in the same device using the first-order correction in 
DRAG, by switching from rotations around a single axis 
induced by Gaussian-modulated microwave tones to rota- 
tions performed using Gaussians and derivatives of Gaus- 
sians applied on two perpendicular axes. We tune up the 
pulses using a set of calibration rotation experiments with 
results that agree with the model outlined in Ref. [15] . 
Randomized benchmarking (RB) is enrployed to show the 
reduction of the average single-qubit gate error [H H51H7] . 
For the shortest gate width of 4 ns, we find an improve- 
ment by a factor of ~ 15 down to a minimum EPG of 
0.007 ± 0.005, which is at the limit imposed by two-level 
relaxation. Independently, the DRAG technique is also 
being implemented at UCSB with superconducting phase 
qubits [18) . 

The optimal control technique of DRAG in Ref. [15] 
prescribes a simple pulse-shaping protocol for reducing 
single-qubit gate errors due to the presence of a third 
level. Neglecting the cavity, which is detuned away from 
any transitions, the driven three-level system, or qutrit, 
is described by the Hamiltonian 
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where <rj = \j — 1) (j\ and g^ = \j) (j — 1| are lower- 
ing and raising operators, w,- j_i are transition energies 
with the ground state energy set to zero, and Ai = 1, 
A2 = A = ill. 2/^0,1.) gi ye the relative drive coupling 
strengths of the o 1 and 1 o 2 transitions, fio.i 
and fii 2j respectively. We can define the anharmonicity 
of the system as the difference between the -f-> 1 and 
1 2 transition frequencies, a± = wj. 2 — Wo,i. The drive 
is only turned on for a fixed gate duration t g and given 
by £{t) = £ x (t) cos(w d i) + £ y {t) sin(cj d i), where £ x > y are 
independent quadrature controls. 

Given a large a\ 3> Wo,i, the effective Rabi drive 
rate J7i 2 to induce any direct or time- dependent tran- 



sitions to |2) can be negligible. However, for a sys- 
tem such as the transmon, an is only ~ 3 — 5% of 
Wo,i- There are two specific qubit gate errors which 
arise due to this reduced anharmonicity. With a non- 
negligible fli,2, it becomes possible for the gate targeting 
the -H- 1 transition to directly populate |2), leaving 
the qubit subspace. However, a second and more dom- 
inant error is a temporary population of |2) during the 
course of a control pulse to the O 1 transition, lead- 
ing to the addition of a phase rotation to the intended 
gate. Although Gaussian control pulses (characterized 
by a width a) are often the paradigm due to their local- 
ized frequency bandwidths given by B = l/2na, leak- 
age errors can occur as gate times are reduced such that 
B is comparable to a.\. A simplified correction proto- 
col to the leakage errors as prescribed in Ref. [T5] is to 
apply an additional control on the quadrature channel, 
£ v (t) — (3£ x (t) and a dynamical detuning of the drive 
frequency S(t) = £ x (t) 2 (— 4/?ai + A 2 )/4ai, where f} is a 
scale parameter. For a qutrit driven without dynamical 
detuning [T!5], the optimal /3 = A 2 /4ai. 

The experiments are performed in a circuit QED sam- 
ple consisting of two transmons coupled to a coplanar 
waveguide resonator. The sample fabrication and experi- 
mental setup are described in Ref. [7] ■ The two transmons 
(designated L and R) are detuned from one another with 
ground to excited state (0 ■<-> 1) transition frequencies of 
u) i /2tt = 8.210, 9.645 GHz and the ground state cavity 
frequency is ujc/2tt = 6.902 GHz. The anharmonicities 
of the transmons are found using two-tone spectroscopy 
measurements 20J to be «[' = 330, 300 MHz and co- 

T R 

herence times are measured to be T x ' = 1.2, 0.9 ^s and 
T 2 * L ' R = 1.5, 1.1 fjs. 

To implement DRAG to first order and perform single- 
qubit gates on the transmons, we use an arbitrary- 
waveform generator to shape microwave-frequency pulses 
with quadrature control, permitting rotations about ei- 
ther the x- or y-axis of each qubit. We fix the drive 
frequency to wo,i, and the pulse amplitudes and phases 
define the rotation angle and axis orientation, respec- 
tively. When performing an x rotation, £ x (t) is a Gaus- 
sian pulse shape [Fig.HTa)], while the derivative of the 
Gaussian is applied simultaneously on the other quadra- 
ture, £ y (t) = j3£ x (t). All of the pulses are truncated to 
2<T from the center with an added buffer time of 5 ns to 
ensure complete separation between concatenated pulses. 

A simple test sequence is used to tune up the scale 
parameter j3 as well as to demonstrate the effect of 
using first-order DRAG pulses versus standard single- 
quadrature Gaussians. The sequence consists of pairs 
of -k and 7r/2 pulses around both the x- and j/-axes. An 
important feature of this sequence is that the final av- 
erage z-projection of the single qubit, (cr z ), will ideally 
take on values from the set S = {+1, 0, —1}, making any 
deviations easily visible. 
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FIG. 1. (color online) (a) Gaussian and derivative pulse 
shapes applied to the in-phase and quadrature control chan- 
nels, respectively, for implementing DRAG to first order. 
(b) Measured (a z ) on qubit L (shaded red bars) with Gaus- 
sians with a — 1 ns for a test sequence consisting of pairs of 
7r and 7r/2 rotations. The slash filled bars correspond to a 
master-equation simulation of a three-level system with pa- 
rameters of the sample tested. The grey filled bars reflect 
ideal values, (c) Similarly measured (a z ) on qubit L but us- 
ing derivative pulses on the quadrature with a scale factor 
/? L = 0.4 (shaded red bars), also overlaid on the ideal values 
(shaded dark grey bars). 



Using Gaussian pulse shaping with a = 1 ns and im- 
plementing the test sequence for qubit L, we find sig- 
nificant deviations from S, as shown in the solid red 
bars of Fig.HTb). The theoretical results for each pair 
of rotations are shown with solid grey bars in the back- 
ground. Note that there are some concatenated rotations 
for which \{cr z )\ > 1> which is because the measurement 
calibration is performed using a 7r pulse which is itself 
non-ideal and plagued by the same errors induced by the 
third level. We observe that the largest errors occur when 
the two rotations are around different axes, which indi- 
cates the presence of significant phase error, or a residual 
z rotation after the first gate. 

We repeat the same test sequence, but applying the 
derivative of the Gaussian to the quadrature channel. By 
varying /3, it is possible to find an optimal value such 
that the measurements of (a z ) agree very well with the 
theoretical predictions. The shaded red bars of Fig.HTc) 
show measured (<r z ) for qubit L using j3 = 0.4. Here, 
deviations from the ideal grey bars decrease to < 2%. We 
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FIG. 2. Randomized benchmarking for qubit L with a — 1 ns 
using (a) Gaussian pulses, and (b) additional Gaussian deriva- 
tive pulses on the quadrature channel. The scattered grey 
points are extracted fidelities for 32 RB sequences, truncated 
at different numbers of gates. A remarkable reduction in the 
extracted average error per gate (black squares) of the bench- 
marking results is observed going from (a) to (b). The error 
bars indicate the variance of all the RB sequences and are 
smaller than the squares in (b). 



have also applied the DRAG protocol for qubit R, finding 
the optimal value /? = 0.25 (data not shown). From the 
experimental determination of /? and ai, we can infer the 
second excited state coupling strengths A L,R = 1.82, 1.41. 
Using A L and the three-level model of Eq.[l] a master 
equation simulation for the Gaussian shaping gives the 
red hash-filled bars in Fig.HTa), which demonstrates good 
agreement with the experiment. 

We find excellent agreement for A with a simple calcu- 
lation for the anticipated A in a cavity-transmon coupled 
system. The cavity modifies the drive strengths Slrj.i and 
^1.2 due to its filtering effect. Specifically, for a transmon 
in a cavity, we have 



%-ij =£ /] 
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nu>c - Wj-i. 



(2) 



where n indexes the cavity mode, j = 1, 2 for the trans- 
mon excitation level, 7 = {1, —1} depending on whether 
the qubit is located at the input or output side of the 
cavity, and gi t j is the matrix element coupling the i ■<->• j 
transmon transition to the cavity [T^] . Using the relevant 
parameters of the two transmons in the experiment and 
including only the fundamental mode of the cavity, we 
find A L ' R ) = 1.85, 1.57, within 12% of those determined 
from the test sequence. There are corrections to the drive 
due to the higher modes of the cavity, but it is difficult to 
use Eq. pi) to estimate as a result of cutoff dependence. 
We characterize the degree of improvement to single- 



qubit gates by using the technique of randomized bench- 
marking (RB) [16] . RB allows us to determine the av- 
erage error per gate through the application of long se- 
quences of alternating Clifford gates (Rx, y ) and Pauli 
gates, chosen from {t,RZ,R*,R%} [21]. We use the RB 
pulse sequences originally given in Ref. [TB] and adapted 
to superconducting qubits in Ref. [4] for both the Gaus- 
sian and the derivative pulse shaping for transmon L. We 
truncate the randomized sequences at various lengths and 
compare the resulting measurement of (a z ) to the ideal 
final state to obtain the fidelity J- . There is an expo- 
nential decrease in J- with an increasing number of gates 
in the randomized sequences. This RB protocol is then 
repeated for various pulse widths, corresponding to dif- 
ferent Gaussian standard deviations, a € {1, 2, 3, 4, 6} ns. 

Using the Gaussian shaping, we find a large reduction 
in fidelity with the shortest pulses, a = Ins [Fig.[2ja)]. 
The scattered grey points give J- for 32 different ran- 
domized sequences applied as a function of the number 
of gates in the sequences. When averaged together, we 
observe a simple decay of J 7 as a function of the number of 
gates (solid black squares). Fitting the data with an ex- 
ponential decay (solid black line), we extract an average 
error per gate, EPG = 1 - T of 0.13 ± 0.02. However, 
when employing the first-order DRAG, we find a dra- 
matic improvement in the gate performance at a = 1 ns 
[Fig.l^b)]. There is a significant reduction in the spread 
of the grey points corresponding to all the different ran- 
domized sequences, and a fit (solid black line) to the 
exponential decay of the average fidelity (solid black 
squares) gives EPG = 0.007 ± 0.005. 

Figure[3] summarizes the improvement to EPG for dif- 
ferent a by using DRAG. The solid squares are the 
EPG found using Gaussians, revealing a minimum of 
0.02 ± 0.007 at a = 3 ns, before considerably increasing 
for shorter pulse lengths. Excellent agreement is found 
with a master equation simulation (dashed line) of the 
gate error for a qutrit system incorporating only deco- 
herence times and coupling strengths measured in inde- 
pendent experiments. Using first-order DRAG, we find 
the solid circles in Fig.[3j which follow a monotonic de- 
crease in EPG with decreasing a. Here again we have 
included a master equation prediction (solid line) of just 
a single qubit with the same parameters, also giving ex- 
cellent agreement with the experimentally determined 
values and demonstrating that DRAG has reduced the 
response of the system to be like that of a single qubit. 

Finally, implementing DRAG on both qubits simul- 
taneously, we can also generate and detect higher qual- 
ity two-qubit states. Performing state tomography to 
obtain the two-qubit density matrix p via joint readout 
[111 122] requires 15 linearly independent measurements, 
corresponding to the application of all combinations of J, 
i?J, Rx , and Ry on the two qubits prior to measure- 
ment. Thus, errors in these analysis rotations in addition 
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FIG. 3. (color online) Comparison of single-qubit gate er- 
rors with and without DRAG. Error per gate for the left 
qubit extracted from randomized benchmarking for different 
gate lengths using both Gaussian pulses (red squares) and 
first-order DRAG pulses (blue squares). Excellent overlap 
with theory for gate error including qubit decoherence (black 
curve) suggests that the DRAG pulses successfully eliminate 
the errors due to the presence of higher levels. Gate errors 
down to 0.007, which are otherwise unattainable with Gaus- 
sian pulses, are reached using DRAG. 



to the state preparation pulses can result in incorrect de- 
termination of p. The two-qubit Pauli set P 11 can be 
used to visualize p for the state |1) L £g> |1) R having used 
Gaussian [Fig.[4|a)] and DRAG [Fig.Qb)] pulse shaping. 
P consists of ensemble averages of the 15 non-trivial com- 
binations of Pauli operators on both qubits. The ideal P 
of the state is characterized by unit magnitude in (ZI), 
(IZ) , and (ZZ) and zero for all other elements. We can 
see that with the standard Gaussian pulse shaping, there 
are substantial (~ 50 — 100% of unity) deviations on ide- 
ally zero elements, whereas with the DRAG pulses, the 
Pauli set bars are very close to their ideal values. 

By implementing a simple approximation to the op- 
timal control pulses for a multi-transmon coupled-cavity 
system, we have reduced gate errors below the 10 -2 level, 
limited by decoherence. The agreement of the various ex- 
periments with and without DRAG pulse shaping with a 
qutrit model reflects that gate errors due to the coupling 
to a higher excited state can be minimized while contin- 
uing to shorten gate time. Moving forward with opti- 
mal control, a tenfold decrease in gate time to approach 
~ 1 ns through improved electronics or a tenfold increase 
in coherence times to ~ 10 ps would place us right at 
the quoted 10 -3 fault-tolerant threshold [15 . Further- 
more, DRAG is extendable to systems of more than two 
multi-level atoms for quantum information processing, 
and has already been employed to enhance single-qubit 
gates in a circuit QED device with four superconducting 
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FIG. 4. (color online) Measured two-qubit Pauli sets for 
preparing the state jl,l) with (a) Gaussian pulses and 



(b) DRAG pulses (/3 L 



0.4, /3 



R 



0.25) applied to the 



quadrature channels of both transmon L and qubit R. The 
ideal Pauli set is shown in grey. 



qubits [25] , 
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